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Water  concentration  in  proton  exchange  membrane  (PEM)  fuel  cells  strongly  influences  performance  and 
durability  which  demands  for  fundamental  understanding  of  water  transport  mechanisms.  The  system 
efficiency  can  be  significantly  improved  with  greater  understanding  of  water  flux  dynamics  through  the 
membrane  and  its  dependence  on  the  internal  conditions  of  the  fuel  cell.  Therefore,  a  two-dimensional, 
non-isothermal,  dynamic  model  of  a  100  W  open  cathode,  self-humidified  PEM  fuel  cell  system  has  been 
developed,  that  is  capable  of  representing  system  specific  control  mechanisms  for  water  and  thermal 
management.  The  model  consists  of  three  coupled  submodels  based  on  energy,  momentum  and  water 
mass  balance  of  the  system.  The  work  is  based  on  experimental  observations  of  the  investigated  fuel  cell 
stack,  for  which  the  crucial  coefficients  for  water  transport,  namely  the  diffusion  and  the  electroosmotic 
drag  (EOD)  coefficient  have  been  determined.  The  diffusivity  of  water  vapor  through  the  MEA  at  30  °C  was 
determined  to  be  3.3  x  10~8  m2  s-1  and  increases  by  3  x  10~10  m2  s-1  °C~1  up  to  50  °C  stack  temperature. 
The  EOD  coefficient  was  measured  as  0.47-0.48  water  molecules  per  proton  at  stack  currents  from  1  to 
3  A.  Validation  of  the  steady  state  and  the  dynamic  model  by  using  experimental  data,  directly  obtained 
from  laboratory  tests,  has  shown  that  the  model  predictions  match  the  experimental  data  well. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

For  the  past  20  years  astonishing  progress  in  terms  of  PEM  fuel 
cell  materials,  component  design,  production,  and  system  power 
density  improvements  have  been  achieved.  However,  there  is  still  a 
lot  to  be  done  in  the  field  of  fuel  cell  system  controls,  which  makes  it 
essential  to  understand  the  different  physical  phenomena  within  a 
fuel  cell  and  how  they  need  to  be  controlled  in  order  to  improve  effi¬ 
ciency,  operating  range  and  durability.  The  hypothesis  is  that  if  the 
water  movement  within  a  PEM  fuel  cell  could  be  controlled  quickly 
to  maintain  optimal  membrane  water  content  and  minimal  liquid 
water,  efficiency  would  be  improved.  As  shown  in  the  experiments 
of  Springer  et  al.  [1],  membrane  proton  conductivity  is  a  strong 
function  of  water  content.  Thus,  the  performance  of  PEM  fuel  cells 
is  sensitive  to  membrane  hydration.  Although  water  is  produced 
during  the  reaction,  the  anode  catalyst  layer  is  often  dehydrated 
because  water  is  dragged  from  the  anode  to  the  cathode  by  pro¬ 
tons  moving  through  the  membrane,  which  is  called  electroosmotic 
drag  (EOD). 

Besides  the  EOD,  the  main  water  transport  mechanism  in  a 
PEM  fuel  cell  is  diffusion  through  the  membrane  due  to  concentra¬ 
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tion  differences  between  anode  and  cathode.  The  third  transport 
mechanism  is  hydraulic  permeation,  which  is  caused  by  pressure 
difference.  EOD  always  transports  water  from  the  anode  to  the 
cathode  whereas  diffusion  can  occur  in  both  directions. 

Water  is  needed  to  maintain  good  proton  conductivity  and 
therefore  has  to  be  kept  in  the  membrane,  however  liquid  water 
on  the  catalyst  reduces  the  active  area,  and  in  the  GDL  it  hinders 
the  reactant  gases  from  diffusing  to  the  catalyst  surface  and  thus 
reduces  performance.  The  goal  is  to  maintain  an  optimal  water  con¬ 
centration  in  the  membrane  electrode  assembly  (MEA)  by  keeping  a 
balance  between  the  two  conflicting  requirements.  Thus,  to  control 
water  transport  within  a  fuel  cell  system  and  thereby  optimize  the 
membrane  hydration  at  any  operation  point,  proper  dynamic  water 
management  strategies  have  to  be  developed.  This  has  recently 
been  analyzed  by  Hussaini  and  Wang  [2]. 

In  order  to  characterize,  understand  and  manipulate  the  water 
transport  mechanisms,  experimental  work  is  needed  as  well  as 
a  mathematical  model  that  describes  the  physical  phenomena 
[3]. 

A  2D  isothermal  model  of  the  MEA  of  a  PEM  fuel  cell  including 
the  influence  of  convection  in  the  gas  flow  channels  was  devel¬ 
oped  by  Gurau  et  al.  [5].  This  model  accounts  for  the  concentration 
variations  along  the  interface  between  the  gas  diffusion  layer  and 
the  catalyst  layer,  which  is  related  to  the  gas  transport  in  the  cou¬ 
pled  domain  of  the  gas  flow  channel  and  the  gas  diffusion  layer. 
However,  fluid  dispersion  in  the  porous  media  is  disregarded. 
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Nomenclature 

Symbols 

A  cross-sectional  area  (m2) 

c  concentration  (mol  m-3) 

Cp  specific  heat  capacity  (J  mol-1  K-1 ) 

D  diffusion  coefficient  (m2  s-1 ) 

€  EOD  coefficient 

r]  efficiency 

<2>h2o  water  mass  flux  (kgs-1  m-2) 

0q  heat  flux  ( W  m-2 ) 

HLhv  lower  heating  value  H2  (1<J  mol-1 ) 

I  current  (A) 

J  molar  flux  (mol  s-1  m-2 ) 

k  thermal  conductivity  (W  m-1  K-1 ) 

k  permeability  (m2) 

M  molar  mass  (kg  mol-1 ) 

m  mass  (kg) 

m  mass  flow  rate  (kg  s-1 ) 

fi  dynamic  viscosity  (Pas) 

nbpp  number  of  bipolar  plates 

nceu  number  of  cells 

ncpp  number  of  channels  per  plate 

P  pressure  (Pa) 

Q.  heat  energy  (W) 

p  density  (kg  m-3) 

T  temperature  (I<) 

t  time  (s) 

v  velocity  (m  s-1 ) 

V stack  stack  voltage  (V) 

W  work  (W) 

Xh2o  humidity  mass  ratio 

Subscripts 
act  active 

an  anode 

ca  cathode 

ch  channel 

cons  consumed 

d  dry 

dp  dew  point 

el  electrical 

gen  generated 

in  inlet 

m  measured 

out  outlet 

sc  short  circuit 

tot  total 


The  amount  of  waste  heat  produced  by  a  PEM  fuel  cell  is  similar 
to  its  electrical  power  output,  depending  on  its  voltage.  Moreover  it 
only  tolerates  a  small  temperature  deviation  from  its  design  point 
for  best  performance,  stability  and  durability.  Therefore  a  three- 
dimensional,  non-isothermal  model  was  developed  by  Ju  et  al.  [6]. 
The  model  accounts  for  various  heat  generation  mechanisms  and 
combines  them  with  the  electrochemical  and  mass  transport  mod¬ 
els.  A  three-dimensional,  non-isothermal,  two-phase  flow  model 
was  developed  by  Wang  and  Wang  [7],  which  was  applied  by  Basu 
et  al.  [8]  in  order  to  study  the  phase-change  phenomena  in  the 
cathode  GDL  of  a  PEM  fuel  cell  and  has  finally  been  extended  to  a 
complete  two-phase  model  for  an  entire  PEM  fuel  cell,  including 
two  phase  flow  in  the  gas  channels,  by  Basu  et  al.  [9]. 

Similar  to  the  model  of  Gurau  et  al.  [5],  but  considering  fluid 
dispersion  in  the  porous  media,  two  modeling  modes  of  a  2D 
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Fig.  1.  ID  and  2D  modeling  domains  for  a  PEM  fuel  cell  [4]. 


isothermal  model  have  already  been  implemented  in  COMSOL 
Multiphysics  by  Shi  and  Wang  [10].  As  defined  in  Fig.  1  the  compu¬ 
tational  domain  of  a  2D  model  can  either  be  a  partial  cross-section 
parallel  (x-z-direction)  or  perpendicular  (y-z-direction)  to  the  gas 
flow  direction  in  the  gas  channel.  They-z-model  serves  for  analysis 
of  fluxes  and  concentrations  in  the  gas  diffusion  and  catalyst  layers 
and  includes  the  effect  of  ribs  or  lands  between  the  channels.  This 
model  is  also  used  for  investigating  fuel  cells  with  interdigitated 
flow  patterns.  If  the  field  of  study  concentrates  on  the  analysis  of 
reactant  gases  and  water  vapor  concentration  along  the  channel, 
the  x-z-model  is  preferred.  The  two  models  can  be  combined  to 
describe  the  overall  behavior  of  PEM  fuel  cells  in  all  directions,  as 
shown  by  Shi  and  Wang  [10]. 

The  water  transport  equations  in  most  of  these  models  are 
based  on  the  ex  situ  water  transport  experiments  of  a  Nafion  117 
membrane  performed  by  Springer  et  al.  [1].  These  experiments 
have  created  a  baseline  for  the  industry  and  show  the  relation¬ 
ship  between  EOD  and  water  diffusion  through  the  membrane  with 
respect  to  membrane  water  activity  and  membrane  temperature. 
Even  though  the  diffusion  and  EOD  data  for  the  membrane  were 
accurate,  direct  application  of  such  data  to  a  real  fuel  cell  may  not 
be  appropriate  due  to  the  fact  that  this  data  was  collected  ex  situ, 
and  cannot  be  considered  constant  because  the  coefficients  vary 
significantly  depending  on  the  membrane  type  and  on  the  operat¬ 
ing  conditions,  such  as  temperature.  In  situ  experiments  of  Husar 
et  al.  [  1 1  ]  showed  that  the  EOD  of  a  Nafion  115  membrane  increases 
significantly  with  temperature  and  current  density  and  that  water 
diffusivity  of  membrane  is  lower  than  that  reported  by  Springer 
et  al.  [1]. 

A  recently  published  review  of  water  balance  in  the  MEA  by 
Dai  et  al.  [3]  states  that  further  work  is  needed  to  better  under¬ 
stand  the  fundamentals  of  water  transport  in  the  MEA,  not  only  to 
improve  performance,  but  also  to  develop  new  materials  for  better 
water  management  and  to  improve  durability.  In  order  to  develop 
and  simulate  dynamic  water  management  strategies  that  match 
the  application  load  requirements  and  the  operating  conditions, 
new  models  need  to  be  based  on  a  broad  understanding  of  water 
transport  in  the  MEA. 

This  article  describes  the  developed  dynamic  thermal  and  water 
distribution  model,  as  well  as  the  performed  experimental  work 
and  the  model  validation  of  an  open  cathode,  self-humidified  PEM 
fuel  cell. 

2.  Specific  stack  characteristics 

This  work  treats  the  modeling  of  the  water  and  heat  transfer  of 
the  commercially  available  1 00  W  PEM  fuel  cell  system  H-l  00  from 
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Horizon  Fuel  Cells.  This  open  cathode  system  is  self-humidified  and 
air  cooled.  It  includes  a  cooling  fan  directly  attached  to  the  fuel  cell 
housing,  which  removes  heat  from  the  stack  by  forced  convection 
and  at  the  same  time  provides  oxygen  to  the  cathode.  The  anode 
system  has  electromagnetic  valves  on  both  the  inlet  and  outlet. 
The  outlet  valve  is  usually  closed  and  the  pressure  is  controlled  by 
a  forward  pressure  regulator.  It  mainly  runs  in  a  dead-ended  mode, 
however  a  periodical  hydrogen  purge  removes  water  and  nitrogen 
that  has  crossed  over  from  the  cathode  that  would  otherwise  hin¬ 
der  the  transport  of  reactant  gas  to  the  catalyst  layer.  A  very  quick 
short  circuit  is  applied  to  the  fuel  cell  to  create  water  and  heat  in 
the  cathode  catalyst  layer.  With  the  manufactures’  controller  the 
interval  of  the  hydrogen  purge  and  the  short  circuit  is  indepen¬ 
dent  of  the  stack  conditions.  This  means  that  even  if  the  fuel  cell 
does  not  require  a  purge  or  a  short  circuit,  the  system  performs  it 
anyway,  which  reduces  efficiency.  In  order  to  increase  efficiency, 
without  reducing  the  robustness  or  operating  range  of  the  system, 
a  broad  understanding  of  the  water  transport  inside  the  fuel  cell  is 
necessary,  which  can  be  studied  by  the  developed  model. 


3.  Model  description 

3A.  Modeling  strategy 

The  developed  model  is  used  to  simulate  and  study  the  effects  of 
the  dynamic  control  mechanisms  for  water  management,  namely 
the  fan,  the  periodical  hydrogen  purge  and  the  short  circuit  to  relate 
them  to  the  fuel  cell  performance.  By  controlling  the  concentration 
of  water  vapor  and  the  additional  creation  of  water  due  to  the  short 
circuit,  membrane  hydration  and  fuel  cell  flooding  can  be  managed. 
As  water  distribution  and  transport  is  dependent  on  temperature 
the  model  has  to  include  not  only  the  mass  balance  but  also  the 
energy  balance  of  the  H-100  fuel  cell  system.  Since  this  model  con¬ 
centrates  on  the  water  transport  within  the  cell,  the  current  density 
distribution  is  not  modeled.  In  order  to  facilitate  the  model,  the  cur¬ 


rent  density  at  the  cathode  catalyst  layer  is  assumed  to  be  constant 
in  the  direction  along  the  flow  channel. 

To  visualize  the  simulation  results,  the  mathematical  model  is 
implemented  in  COMSOL  Multiphysics,  a  finite  element  analysis 
software  environment.  By  means  of  COMSOL  the  model  can  be 
solved  numerically  for  the  specific  geometry. 
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Fig.  3.  Model  geometric  subdomains. 
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The  trade  off  between  accuracy  and  computing  capacity  should 
be  considered  in  choosing  the  model  dimension.  Within  this  work 
a  2D  model  is  developed  that  provides  better  simulation  accuracy 
than  ID  models  and  has  a  higher  computational  efficiency  than  3D 
models,  which  also  already  exist  in  the  literature  [6,7,9]. 

Referring  to  Fig.  2,  the  2D  model  describes  water  and  heat  prop¬ 
agation  and  distribution  in  the  x-z-cross-section  of  a  single  cell 
within  the  stack.  The  cross-section  can  either  go  through  a  rib  of  a 
bipolar  plate  or  through  a  channel.  Since  the  primary  objective  of 
the  model  is  to  describe  the  effects  of  water  transport,  the  air  mass 
flow  through  the  channel  plays  an  important  role  and  so  the  sec¬ 
ond  option  for  the  cross-section  location  is  chosen.  A  2D  thermal 
analysis  in  the  y-z-cross-section  has  shown  that  the  temperature 
difference  at  the  channel  walls  in  the  same  plane  is  negligible,  due 
to  the  relatively  high  thermal  conductivity  of  the  solid  sections 
and  the  fully  developed  flow  through  the  channel.  Thus,  the  heat 
removal  through  the  two  land  sides  of  one  channel,  that  do  not 
appear  in  the  2D  model  in  the  x-z-cross-section,  is  assumed  to  be 
equal  to  the  heat  removal  through  the  other  two  channel  walls.  This 
simplification  is  valid  in  this  model,  because  the  channel  is  square 
and  the  thickness  of  the  land  is  equal  to  the  side  of  the  square,  as 
shown  in  Fig.  2.  Considering  also  that  the  x-z-model  is  infinite  in  the 
y-direction,  which  means  that  there  is  an  infinite  channel  without 
any  land,  the  effective  heat  removal  area  in  both  models  is  equal. 
Thus  the  overall  heat  removal  from  the  cell  can  be  modeled  by  just 
one  2D  model  in  the  x-z-cross-section. 

Fig.  3  shows  the  geometry  in  the  x-z-plane  of  the  modeled 
cell  within  the  H-100  stack  referred  to  the  cross-section  marked 
in  red  in  Fig.  2.  The  modeled  sub-domains  for  the  simulation  of 
water  propagation  and  distribution  are  the  cathode  flow  channel, 
the  cathode  GDL  (including  a  microporous  layer)  and  two  anode 
GDL’s  (including  a  microporous  layer).  To  simulate  the  heat  trans¬ 
fer  within  the  cell,  the  Grafoil  gasket  and  the  bipolar  plate  also  have 
to  be  considered,  due  to  conductive  heat  transfer.  The  right  bound¬ 
ary  of  the  bipolar  plate  is  equal  to  the  left  boundary  of  the  cathode 
flow  channel  in  terms  of  heat  transfer,  because  only  one  repeating 
unit  of  the  stack  is  modeled. 

In  terms  of  modeling,  the  material  properties  and  dimensions  of 
the  components,  as  well  as  the  anode  channel  configuration  have 
to  be  known.  The  thicknesses  of  the  different  layers  are  listed  in 
Table  1.  The  channel  length  is  25  mm. 

Since  the  material  of  the  membrane  is  unknown,  a  bulk  diffu¬ 
sion  coefficient  for  the  whole  MEA  was  determined  experimentally. 
The  MEA  on  the  cathode  side  includes  the  GDL  and  the  microp¬ 
orous  diffusion  layer  and  the  cathode  catalyst  layer.  On  the  anode 
side  the  MEA  includes  the  membrane,  anode  catalyst  layer,  micro¬ 
porous  diffusion  layer,  and  two  GDL’s.  To  clarify,  the  membrane  is 
included  in  the  first  anode  GDL,  which  is  denoted  by  the  dashed 
line  in  Fig.  3.  Thus,  in  terms  of  water  transport  and  generation  the 
more  important  cathode  catalyst  boundary  still  remains. 

According  to  the  model  geometry,  Fig.  4  shows  a  schematic  of 
the  different  physical  phenomena  that  occur  within  the  fuel  cell, 
how  they  are  coupled  to  each  other  and  how  they  are  treated  in  the 
model. 


Table  1 

H-100  component  thicknesses. 


Component 

Thickness  (mm) 

Bipolar  plate  web 

0.70 

Grafoil  gasket 

0.55 

Anode  GDL  (flow  channel) 

0.20 

Anode  GDL  (cover) 

0.20 

Microporous  layer 

0.20 

Membrane 

0.05 

Cathode  GDL 

0.20 

Cathode  flow  channel 

1.50 

The  heart  of  the  model  is  the  water  transport  submodel,  which 
describes  the  distribution  of  water  vapor  concentration  in  the  MEA 
and  the  flow  channels.  In  order  to  describe  convective  transport,  the 
water  transport  submodel  is  coupled  to  the  momentum  transport 
sub  models  of  anode  and  cathode,  which  themselves  are  linked  with 
the  water  transport  submodel  because  the  density  of  the  reactant 
gases  is  dependent  on  the  amount  of  water  in  the  gas.  Since  density 
is  also  dependent  on  temperature,  the  energy  transport  submodel 
has  to  be  included,  as  well.  As  shown  in  Section  4,  the  diffusion  coef¬ 
ficient  is  also  a  function  of  temperature,  which  links  the  diffusive 
mass  transport  to  the  temperature  distribution.  The  water  trans¬ 
port  from  anode  to  cathode  due  to  the  electroosmotic  drag  effect  is 
a  function  of  the  stack  current,  which  is  set  by  the  external  load. 

The  energy  transport  model,  which  describes  the  temperature 
distribution  within  the  cell  includes  not  only  conductive  heat  trans¬ 
fer  through  the  MEA,  the  Grafoil  gasket  and  the  bipolar  plate,  but 
also  the  convective  cooling  by  the  fan.  Therefore  it  is  coupled  to  the 
velocity  field  in  the  cathode  channel,  obtained  by  the  momentum 
transport  submodel,  as  described  in  Section  3.2.2. 

Depending  on  the  external  load,  a  certain  amount  of  heat  and 
water  is  generated  at  the  cathode  catalyst  surface  due  to  the  elec¬ 
trochemical  reaction.  This  generation  can  be  described  by  the  stack 
current  and  voltage  that  is  set  by  the  external  load.  The  resulting 
fluxes  of  water  and  heat  are  treated  as  an  input  to  the  water  and 
energy  transport  submodel,  respectively.  The  different  submodels 
are  explained  in  the  following  sections. 

A  further  submodel  that  describes  charge  transfer  and  polariza¬ 
tion  curve  affected  by  the  different  types  of  voltage  losses  is  not 
included  in  this  work,  because  the  needed  parameters  to  deter¬ 
mine  voltage  losses,  such  as  exchange  current  density  and  charge 
transfer  coefficient,  are  difficult  to  validate  with  in  situ  testing. 

This  model  considers  that  water  enters  and  exits  the  fuel  cell  in 
the  vapor  state  on  both  anode  and  cathode.  Regarding  the  cathode, 
water  vapor  that  enters  the  cell  has  the  relative  humidity  of  the 
environment.  On  the  cathode  outlet  there  is  always  a  stoichiom¬ 
etry  greater  than  ten  that  does  not  allow  the  gas  to  reach  50% 
relative  humidity  at  an  outlet  gas  temperature  of  approximately 
45  °C,  which  has  also  been  verified  by  experiments.  Thus,  a  single¬ 
phase  model  on  the  cathode  side  is  accurate.  The  only  controversial 
part  would  be  the  anode  where  the  water  could  condense.  Since  the 
anode  reaction  is  orders  of  magnitude  faster  than  the  cathode  reac¬ 
tion  and  hydrogen’s  diffusivity  is  much  higher  than  oxygen’s,  the 
effects  of  liquid  water  on  the  anode  can  be  neglected  for  the  sake 
of  model  simplification.  So  the  model  properly  predicts  water  and 
temperature  distribution  in  the  flow  channels  and  the  GDL,  even 
without  including  two-phase  flow. 

Table  2  shows  how  the  COMSOL  application  modes  are  applied 
to  the  specific  subdomains.  These  application  modes  are  linked 
to  the  different  submodels,  which  are  explained  in  the  following 
section. 


3.2.  Mathematical  model  description 
3.2 A.  Energy  transport 

The  amount  of  energy  brought  into  the  fuel  cell  system  is  given 
by  the  lower  heating  value  of  hydrogen,  because  it  is  assumed  that 
all  of  the  product  water  leaves  the  stack  as  vapor.  The  energy  output 
is  split  into  electrical  energy  and  heat,  thus  the  energy  balance  of 
the  system  is  given  by: 

Wh2  =  Wel  +  Qtot  (1) 

The  total  generated  heat  Qtot  can  be  determined  using  the  fuel 
cell  efficiency,  which  is  defined  by  the  ratio  of  the  energy  output  and 
the  energy  input  of  the  system.  According  to  Eq.  (1 )  this  is  the  ratio 
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Fig.  4.  Model  schematic. 


of  useful  electrical  work  and  the  work  of  the  consumed  hydrogen:  generated  heat  energy  by  including  Eq.  (4)  into  (1 ): 


(2) 


The  work  of  the  consumed  dry  hydrogen  at  25  °C  can  be  calcu¬ 
lated  using  hydrogen’s  lower  heating  value  [4]: 


Wh2  = 


Hlhv 
2  F 


'  f  '  ftcell 


[W] 


(3) 


Combining  Eqs.  (2),  and  (3)  leads  to  the  fuel  cell  stack  efficiency: 


“Hun, stack  = 


Yj. stack 

1 .254  •  ncen 


(4) 


Current  and  voltage  values  that  are  used  to  calculate  efficiency 
are  obtained  by  in-house  experimental  data  of  the  studied  fuel  cell 
stack.  The  stack  efficiency  can  then  be  used  to  describe  the  total 


dot  =  Wh2  -  Wel  =  We,  ■  (  — 2 - l)  [W]  (5) 

y  VLHV,  stack  J 

Assuming  that  all  the  generated  heat  energy  is  released  on  the 
cathode  catalyst  layer,  the  heat  flux  through  this  boundary  is  deter¬ 
mined  by  the  generated  heat  divided  by  the  active  area: 

[Wm-2]  (6) 

Heat  is  transferred  by  two  mechanisms  within  the  fuel  cell.  The 
cathode  flow  channel  is  dominated  by  forced  convection  due  to 
the  fan,  whereas  the  heat  transfer  through  the  MEA  and  the  bipolar 
plate  is  mainly  by  conduction,  which  can  be  seen  in  the  heat  transfer 
simulation  results  in  Fig.  5.  Although  there  is  also  convection  due  to 


Table  2 

COMSOL  application  modes  used  in  the  different  model  subdomains. 


General  heat  transfer 

Incompressible  Navier-Stokes 

Darcy’s  law 

Convection  and  diffusion 
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X 

X 
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X 
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Fig.  5.  Typical  temperature  distribution  in  the  cell  and  the  cathode  flow  channel  ( x  and  y  axis  dimensions  are  in  meters). 


the  diffusing  reactants,  this  transport  can  be  considered  negligible 
because  of  the  high  thermal  conductivity  of  the  GDL  material  [6]. 

Several  thermal  analyses  and  models  of  non-insulated  PEM  fuel 
cells  e.g.  by  Mller  and  Stefanopoulou  [12]  have  shown  that  the  heat 
loss  to  the  environment  by  convection  and  radiation  is  between 
10%  and  20%  of  the  total  waste  heat.  Since  the  fuel  cell  stack  used 
to  validate  the  model  is  encased  in  a  plastic  support  structure,  it 
is  assumed  that  the  heat  loss  to  the  environment  is  minimal  and 
thus  neglected.  Therefore,  only  heat  removal  from  the  system  by 
convective  heat  transport  through  the  cathode  gas  channel  is  con¬ 
sidered. 

Using  conservation  of  energy,  the  heat  equation  that  describes 
conduction  and  convection  within  the  fuel  cell  results  in  [4]: 

(pCp)  •  J  +  (pCp)  •  v  ■  VT  =  V(k .  VT)  +  Se  (7) 

The  first  term  of  the  left-hand  side  describes  the  heating  rate  and 
the  second  term  the  heat  flux  due  to  convection.  The  first  term  on 
the  right-hand  side  accounts  for  conductive  heat  flux  through  the 
media  with  the  thermal  conductivity  k,  described  by  Fourier’s  law. 
The  source  term  Se  represents  heat  flux  due  to  other  heat  sources  or 
sinks,  which  in  this  case  is  the  source  heat  flux  through  the  catalyst 
surface,  described  by  Eq.  (6). 

The  velocity  vector  v  is  obtained  by  solving  the  Navier-Stokes 
equation  for  momentum  transport  within  the  cathode  flow  chan¬ 
nel,  which  is  explained  in  Section  3.2.2. 

The  input  parameters  of  the  model  are  the  inlet  air  temperature 
of  the  cathode,  which  is  necessary  to  determine  the  heat  removal 
from  the  cell  by  convection,  as  well  as  the  cell  voltage  and  current 
that  are  used  to  calculate  heat  generation  at  the  cathode  catalyst  as 
already  described  in  Section  3.2.1 .  The  submodel  schematic,  includ¬ 
ing  all  inputs  and  outputs,  is  shown  in  Fig.  6. 

The  energy  model  is  coupled  with  every  other  sub-model 
because  temperature  has  a  significant  effect  on  the  fluid  properties, 
such  as  densities  of  the  reactant  gases,  their  dynamic  viscosities  and 
also  on  the  diffusivity  of  water  through  the  MEA. 

In  the  H-100  system  heat  is  produced  due  to  the  electrochem¬ 
ical  reaction,  as  well  as  periodically  due  to  the  short  circuit  that  is 
applied  every  10  s,  as  described  in  Section  2. 

During  a  short  circuit  the  useful  electrical  energy  is  zero.  This 
means  that  all  the  energy  is  transformed  into  heat,  which  increases 
temperature  inside  the  fuel  cell  stack.  Nevertheless  some  energy  is 


lost  due  to  circuit  and  contact  resistance,  but  is  neglected  at  this 
point.  Thus,  knowing  the  short  circuit  current,  as  defined  later  in 
Section  3.2.4,  the  released  heat  of  one  cell  can  be  calculated  using 
Eq-(3): 

wh2  =  ■  IsM  =  Qsc  [W]  (8) 

3.2.2.  Cathode  momentum  transport 

In  order  to  describe  convective  mass  transport  and  heat  transfer 
through  the  cathode  flow  channel  forced  by  the  fan,  the  momentum 
transport  has  to  be  determined. 

The  Reynolds  number  for  the  H-100  cathode  flow  channels  at 
maximum  flow  velocity  turned  out  to  be  120.  Knowing  that  the  air 
flow  through  the  channels  is  laminar  and  that  the  pressure  differ¬ 
ence  along  the  channels  and  the  change  in  fluid  density  are  also 
very  small,  the  Navier-Stokes  equation  for  an  incompressible  fluid 
can  be  used  to  model  the  momentum  transfer  through  the  cathode 
flow  channels  [13]: 

^  (If  +*'■ vi')  =  v{-w+/4v?+(v^)t])  (9) 

By  solving  this  equation  numerically,  the  velocity  field  in  the  cath¬ 
ode  flow  channel  can  be  obtained.  Fig.  7  shows  a  schematic  with  all 
input  and  output  parameters  of  the  submodel. 

Since  the  velocity  is  dependent  on  the  gas  density,  which  is  a 
function  of  temperature,  pressure  and  water  concentration,  this 
model  has  to  be  coupled  with  the  energy  and  water  transport  sub¬ 
model.  The  boundary  conditions  for  the  model  are  an  inlet  velocity 
Vo  and  a  constant  outlet  pressure. 

As  the  cathode  air  flow  determines  the  amount  of  heat  that  is 
removed  from  the  stack,  it  has  to  be  controlled  by  the  fan  accord¬ 
ing  to  the  stack  temperature.  Control  strategies  can  be  tested  in 
the  model  simply  by  changing  the  inlet  boundary  condition  from  a 
constant  velocity  to  the  desired  temperature  dependent  boundary 
expression,  which  is  then  coupled  to  the  energy  transport  model. 

3.2.3.  Anode  momentum  transport 

As  already  described  in  Section  2  the  anode  is  purged  periodi¬ 
cally  every  10  s.  This  results  in  a  convective  flux  through  the  anode 
GDL  which  removes  product  water  from  the  GDL  and  the  cata¬ 
lyst  layer.  Because  this  pressure  driven  convective  flux  occurs  in 


S.  Strahl  et  al.  /  Journal  of  Power  Sources  196(2011)  4251-4263 


4257 


Fig.  6.  Schematic  of  the  heat  transfer  model. 


a  porous  media,  the  velocity  field  can  be  calculated  by  Darcy’s  law 
[13]: 


where  k  denotes  the  hydraulic  permeability  of  the  porous  medium 
and  ft  the  dynamic  viscosity  of  the  fluid. 

The  anode  inlet  gas  velocity  is  determined  by  the  measured  mass 
flow  rate.  Thus,  knowing  the  velocity  and  pressure  difference  from 
inlet  to  outlet,  the  permeability  of  the  GDL,  including  the  microp- 
orous  layer,  turned  out  to  be  10-12  m2,  which  accords  to  the  work 
of  Shi  et  al.  [  14]  and  also  is  used  in  the  models  of  Shi  and  Wang  [15] 
and  Meng  and  Wang  [16].  The  calculated  velocity  vector  is  coupled 
to  the  water  transport  model  in  the  anode  GDL  during  a  purge  to 
describe  the  convective  mass  transport.  Fig.  8  summarizes  the  inlet 
and  outlet  parameters  of  the  model. 

Just  as  in  the  cathode  flow  model,  the  anode  flow  model  also  has 
to  be  coupled  to  the  energy  and  water  transport  model  due  to  the 
changes  in  gas  density. 

3.2.4.  Water  transport 

In  order  to  determine  the  water  transfer  rate,  which  is  required 
for  proper  water  management  in  the  fuel  cell  and  to  validate  exper¬ 
imental  and  modeling  work,  the  water  mass  balance  across  the  fuel 
cell  is  needed.  The  full  water  mass  balance  equation  is: 

^H20, ca,in  T  IT1r20 ,an,in  T  rilu20,gen  —  f^H20, ca, out  T  l3lR20,an,out  (11) 

This  equates  the  water  that  enters  and  is  generated  in  the  fuel 
cell  to  the  water  that  leaves  the  fuel  cell.  Water  that  enters  and  exits 
the  cell  at  the  cathode  is  assumed  to  be  in  the  vapor  form,  as  already 
explained  in  Section  3.1 .  The  only  controversial  mass  flow  would  be 


the  anode  outlet,  where  the  water  could  condense.  However,  this 
is  solved  by  placing  a  gas  line  heater  at  the  exit  of  the  anode,  which 
heats  the  gas  up  to  about  70  °C  before  measuring  the  dew  point 
temperature  which  allows  for  the  measurement  of  all  the  water 
leaving  the  anode.  Consequently,  the  mass  flow  of  water  is  in  the 
vapor  form  at  the  location  where  the  dew  point  temperatures  are 
measured. 

Generally  the  mass  flow  rate  of  water  vapor  can  be  expressed  as 
a  fraction  of  the  dry  gas  mass  flow  rate: 

mH2o,ij  =  mUjd  •  xH2o,,j  [kg S_1 1  (12) 

where  i  stands  for  anode  or  cathode  and  j  for  inlet  or  outlet.  The 
determination  of  the  five  different  terms  in  Eq.  (11),  according  to 
the  measurable  variables  of  the  test  station,  is  explained  below. 

1 .  Anode  inlet:  The  total  anode  inlet  mass  flow  rate  is  a  sum  of  the 
hydrogen  mass  flow  rate  and  the  water  vapor  mass  flow  rate,  if 
hydrogen  is  humidified: 

™an,in  —  ^R2,an,in  +  ^R20, an, in  [kgS  ^]  (13) 

Knowing  the  amount  of  water  vapor  entering  with  the  hydrogen 
by  measuring  the  dew  point  temperature  and  the  inlet  mass  flow 
rate  of  dry  hydrogen,  the  anode  inlet  mass  flow  rate  of  water 
vapor  can  be  calculated  by  Eq.  (12): 

^H20, an, in  —  2,an,in  '  *H20, an, in  [kgS  ]  (14) 

2.  Cathode  inlet:  The  total  cathode  inlet  mass  flow  rate  is  a  sum  of 
oxygen,  nitrogen  and  water  vapor  mass  flow  rate: 

™ca,in  —  ™02,ca,in  +  2,ca,in  +  ™R20,ca,in  [kgS  ]  (15) 


Fig.  7.  Schematic  of  the  cathode  flow  model. 


4258 


S.  Strahl  et  al.  /  Journal  of  Power  Sources  196(2011)  4251-4263 


Fig.  8.  Schematic  of  the  anode  flow  model. 


Moreover  a  mass  flow  rate  of  a  gas  stream  can  be  described  by 
the  gas  density,  the  velocity  and  the  cross-sectional  area,  through 
which  gas  flows.  As  the  cathode  inlet  velocity  is  measured  within 
the  test  station  the  cathode  inlet  mass  flow  rate  is  described  by: 

™ca,in  —  Vca,in  '  ^ca  '  Pea, in  [kgS  ]  (16) 

where  Aca  is  the  cross-sectional  area  of  the  fuel  cell  housing 
structure  within  the  environmental  chamber.  Assuming  that  the 
median  flow  velocity  vca  in  is  constant  for  wet  and  dry  air,  because 
the  change  in  mass  flow  rate  results  from  the  change  in  density 
by  adding  water,  the  cathode  inlet  mass  flow  rate  of  water  vapor 
can  then  be  calculated  by  combining  Eqs.  (16)  and  (12): 

™H20,ca,in  —  ™-air,in,d  '  *H20 ,ca,in  —  P ca,in  '  ^ca  •  Pair,d  '  *H20 ,ca,in 

[kgs"1]  (17) 


As  the  cathode  inlet  air  velocity  is  measured  in  a  housing  struc¬ 
ture  outside  the  fuel  cell  stack,  but  has  to  be  used  to  model  a 
single  channel,  the  velocity  has  to  be  adapted  to  the  smaller  flow 
channel  dimensions: 


“ housing 


nbpp  '  ncpp  •  Ach 


■Pm 


[ms 


-ii 


(18) 


3.  Anode  outlet:  The  total  anode  outlet  mass  flow  is  defined  as: 


Wlan,out  —  thn  2, an, out  +  ftlH20 , an, out  [kgs  ’]  (19) 

This  mass  flow  rate  of  water  vapor  can  be  calculated  similar  to 
the  inlet  flow  rate,  but  adding  the  fact  that  hydrogen  is  consumed 
inside  the  fuel  cell.  Therefore  the  anode  outlet  mass  flow  rate  of 
hydrogen  is: 

lPlH2,an,out  —  i^-H2,an,in  ~  ^H2, an, cons)  [kgS  M  (20) 

Using  Faraday’s  law,  the  consumed  hydrogen  mass  flow  rate 
of  one  cell  is  found  by: 

™H2,an,cons  =  ^  [kgS  M  (21) 

Combining  Eqs.  (12)  and  (20)  the  anode  outlet  mass  flow  rate 
of  water  vapor  is  described  by: 

ft*H20,  an,  out  =  (mH2  ,an,in  ^H2, an, cons)  '  *H20, an, out  [kgs  '] 

(22) 


4.  Cathode  outlet:  Regarding  that  only  oxygen  is  consumed  at  the 
cathode  and  the  amount  of  nitrogen  stays  the  same,  the  total 
cathode  outlet  mass  flow  rate  results  in: 


1TT-ca,out  —  (ftlo 2,ca,in  ~  ™02,ca,cons)  +  ™-N2,ca,in  +  ftlH20 ,ca,out 

[kgs-1]  (23) 


The  consumed  oxygen  mass  flow  rate  of  one  cell  can  be 
calculated  using  Faraday’s  law  similar  to  the  consumption  of 
hydrogen.  The  only  difference  is  that  the  number  of  electrons 
per  molecule  of  O2  is  4. 

mo2,ca,cons  =  '  [kg  ]  (24) 

According  to  Eq.  (12),  the  cathode  outlet  mass  flow  rate  of 
water  vapor  is  determined  by  the  following  equation: 

Wh20, ca, out  =  W air, out, d  '  XH20,ca, out  [kgS  ^  ]  (25) 

where  the  outlet  mass  flow  rate  of  dry  air  can  be  determined 
by  subtracting  the  mass  flow  rate  of  consumed  oxygen  from  the 
inlet  mass  flow  rate  of  dry  air: 

P^ai r,out,d  —  ™-air,in,d  ~  ^O 2,ca,cons  [kgS  ]  (26) 

5.  Generated  water:  The  water  that  is  generated  in  the  fuel  cell  is 
a  product  of  the  reaction  of  hydrogen  and  oxygen  and  is  directly 
proportional  to  the  current  passing  through  the  fuel  cell  stack, 
based  on  Faraday’s  law.  Thus,  the  mass  flow  rate  of  generated 
water  results  in: 

Mh2o  r  M  _i, 

m\\2o,gen  =  ~^r  1  [kgs1]  (27) 

As  water  is  only  produced  at  the  cathode  catalyst  surface,  the 
cathode  catalyst  boundary  condition  in  the  2D  model  is  given  by 
the  flux  of  generated  water  through  the  active  area: 

0H2O,gen  =  [kgs”1  ITT2]  (28) 

Table  3  depicts  the  different  contributions  to  the  water  mass 
balance  Eq.  (1 1 )  at  3  A.  As  it  can  be  seen,  the  majority  of  the  water 
enters  and  leaves  the  stack  at  the  cathode  side,  which  is  relatively 
large  compared  to  the  generated  water  at  this  current. 

This  water  mass  balance  is  used  to  determine  the  bulk  water 
vapor  diffusion  coefficient  Dh2o  of  the  MEA,  as  described  in  Section 
4  and  to  validate  experimental  results. 

Since  water  in  a  fuel  cell  is  transported  by  convection,  diffusion 
and  is  also  generated  within  the  cell,  the  different  equations  for 


Table  3 

Contributions  to  the  water  mass  balance  equation  at  3  A. 


Mass  flow  rate 

Absolute  value  (mg  s-1 ) 

mH20,  an,  in 

0.40 

rilH20,ca,in 

60.12 

hlH20 , an, out 

-0.52 

AlH20  ,ca,out 

-65.60 

jtlH20  gen 

5.60 
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Fig.  9.  Typical  water  concentration  distribution  for  constant  anode  and  cathode  flows  at  5  A  (x  and  y  axis  dimensions  are  in  meters). 


each  transport  mechanisms  have  been  combined  in  a  mass  balance 
equation  [13]: 


-  +  V(-Dh2oVch2o)  =  V/h2o  -  V  ■  Vch2o  (29) 


The  first  term  on  the  left-hand  side  of  the  equation  corresponds 
to  the  accumulation  of  water  in  the  system.  The  second  term 
accounts  for  the  diffusive  transport  within  the  MEA,  described  by 
Fields  law.  The  first  term  on  the  right-hand  side  represents  a  source 
flux  of  water  due  to  the  chemical  reaction  and  also  the  EOD.  Finally, 
the  second  term  on  the  right-hand  side  accounts  for  the  convec¬ 
tive  transport  due  to  a  velocity  field  v.  In  the  cathode  flow  channel 
this  field  is  obtained  by  coupling  the  Navier-Stokes  momentum 
transport  to  the  equation  system,  as  described  in  Section  3.2.2, 
whereas  in  the  anode  GDL  the  momentum  transport  during  a  purge 
is  described  by  Darcy’s  law,  as  explained  in  Section  3.2.3. 

Water  is  also  transported  by  the  electroosmotic  drag  from  anode 
to  cathode.  This  flow  rate  can  also  be  described  by  Faraday’s  law, 


because  the  EOD  is  proportional  to  current  [4]: 

^h2o,eod  =  £  •  — y~  •  f  [kgs-1]  (30) 

The  EOD  coefficient  e  represents  the  number  of  water  molecules 
that  are  dragged  from  anode  to  cathode  per  proton.  This  coefficient 
is  determined  experimentally,  as  described  in  Section  4.  Regard¬ 
ing  the  model,  the  EOD  is  treated  as  an  internal  source  flux  at  the 
cathode  and  a  sink  flux  at  the  anode.  A  typical  water  concentration 
distribution  for  constant  cathode  and  anode  flows  is  shown  in  Fig.  9, 
where  the  streamlines  describe  the  flux  of  water  from  the  cathode 
catalyst  layer  to  the  cathode  and  anode  flow  channel.  As  the  H-100 
is  an  open  cathode  stack,  the  cathode  inlet  conditions  are  ambient. 
At  the  anode  dry,  pure  hydrogen  enters  without  humidification. 

Fig.  10  summarizes  all  input  and  output  parameters  of  the 
model,  as  well  as  the  influences  of  the  other  submodels. 

During  a  short  circuit  heat  and  water  are  produced  at  the  cath¬ 
ode  catalyst  layer.  This  mass  flow  rate  of  generated  water  during  a 
short  circuit  can  be  calculated  using  Eq.  (27)  and  replacing  the  stack 
current  by  the  short  circuit  current.  As  the  short  circuit  current  is 


Fig.  10.  Schematic  of  the  water  transport  model. 
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Fig.  11.  I-t- plot  of  a  short  circuit,  (a)  Short  circuit  duration  and  (b)  rising  edge. 


a  function  of  time,  the  root  mean  square  (RMS)  current  is  used  to 
calculate  the  mass  flow  rate 

Mh2o  r  n  _tl 

mH20  ,gen,sc  —  ~2p~  sc  ^  ]  (31) 

where 

Isc  =  VxlOfdt  [A]  (32) 

The  RMS  short  circuit  current  of  the  measured  J-t-curve,  shown 
in  Fig.  11,  results  in  23  A. 

Table  4  lists  the  additional  model  input  parameters,  such  as 
physical  properties,  dimensions  and  constants,  which  are  used  in 
the  model. 

4.  Experimental  determination  of  coefficients 

The  diffusion  coefficient  of  water  vapor  through  the  MEA  (as 
stated  in  Section  3.1  also  includes  all  the  diffusion  layers)  is  depen¬ 
dent  on  temperature  and  water  content.  In  order  to  develop  a 
mathematical  relation  between  diffusion,  temperature  and  water 
content,  water  diffusion  has  to  be  separated  from  the  other  water 
transport  mechanisms,  namely  EOD  and  hydraulic  permeation.  As 
the  experiments  of  Husar  et  al.  [  1 1  ]  have  shown,  water  transfer  due 
to  hydraulic  permeation  is  at  least  an  order  of  magnitude  lower  than 
that  due  to  the  two  other  transport  mechanisms,  and  therefore  can 
be  neglected.  To  separate  diffusion  from  the  EOD,  the  fuel  cell  is 


Table  4 

Physical  parameters  and  properties. 


Description 

Value 

Unit 

Active  catalyst  area,  Aact 

0.00225 

m2 

Bipolar  plate  density,  pppp 

1850 

kg  m-3 

Bipolar  plate  thermal  conductivity,  kppp 

14 

WK-1  m1 

Cross-section  channel,  Ach 

1.5  x  10“6 

m2 

Faraday  constant,  F 

96487 

CmoL1 

GDL  density,  Pgdl 

440 

kg  rrr3 

GDI  permeability,  kGdl 

1  x  10-12 

m2 

GDI  porosity,  eGDL 

0.78 

GDL  thermal  conductivity  inplane,  kGDi,in 

21 

WK-1  m1 

GDL  thermal  conductivity  through-plane, 

k  GDI, through 

1.7 

WK-1  m1 

GDL  thickness,  zGdl 

4x  10-4 

m 

Number  of  bipolar  plates,  ntpp 

21 

Number  of  cells,  nceu 

20 

Number  of  channels  per  plate,  ncpp 

51 

Universal  gas  constant,  R 

8.314 

J  mol-1  K"1 

disconnected  from  the  external  circuit  and  nitrogen  is  used  instead 
of  hydrogen,  which  also  guarantees  that  no  water  can  be  gener¬ 
ated  due  to  possible  crossover  of  hydrogen,  which  would  generate 
water  at  the  cathode  side.  The  experiment  was  carried  out  with  a 
dry  anode  and  a  wet  cathode.  The  diffusive  water  mass  flux  across 
the  membrane  is  the  outlet  water  mass  flux  at  the  dry  side  which 
should  equal  the  difference  between  inlet  to  outlet  water  mass  flux 
of  the  wet  side.  Since  the  H-100  is  an  open  cathode  fuel  cell  it  is 
difficult  to  perform  a  similar  diffusion  experiment  in  the  opposite 
direction,  with  a  wet  anode  and  a  dry  cathode.  However,  random 
tests  at  stable  points  where  the  anode  humidity  was  higher  than  at 
the  cathode  side  have  shown  a  direction  independence  of  the  dif¬ 
fusion  coefficient.  The  dew  point  temperature  of  the  wet  cathode 
was  kept  constant  at  20  °C,  which  means  that  the  partial  pressure 
of  water  vapor  in  air  does  not  change  with  temperature,  unlike 
with  using  a  constant  relative  humidity.  The  ambient  temperature 
was  increased  from  30  to  50  °C  with  a  step  size  of  10°C  by  using 
an  environmental  chamber.  At  each  point  enough  time  was  given 
for  the  anode  dew  point  temperature  to  become  stable,  so  that  a 
steady  state  can  be  assumed.  The  chamber  temperature  then  equals 
the  stack  temperature.  Moreover,  the  maximum  cathode  flow  rate 
was  applied  in  order  to  minimize  the  concentration  difference  from 
inlet  to  outlet  at  the  wet  side. 

Fig.  12  shows  the  membrane  diffusivity  of  water  vapor  at  dif¬ 
ferent  temperatures.  Compared  to  the  work  of  Springer  et  al.  [1] 
the  experimentally  obtained  diffusion  coefficient  is  smaller.  This 
might  be  due  to  the  different  membrane  types  and  thickness  used 


Fig.  12.  Membrane  diffusivity  as  a  function  of  stack  temperature. 
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Fig.  13.  Anode  and  cathode  outlet  water  concentration  as  a  function  of  current. 


Fig.  14.  Dynamic  validation— humidity  fan  test. 


in  this  work.  However,  similar  results  were  found  by  McKay  and 
Stefanopoulou  [17],  who  also  performed  in  situ  measurements  of 
the  diffusion  coefficient  with  a  fuel  cell  stack. 

In  order  to  determine  the  EOD  coefficient  experimentally,  diffu¬ 
sive  water  transport  through  the  membrane  has  to  be  minimized. 
This  can  be  obtained  by  keeping  the  water  concentration  on  both 
sides  equal  and  as  close  to  1 00%  relative  humidity  as  possible,  which 
is  achieved  by  setting  the  same  dew  point  temperature  for  the 
anode  and  cathode  inlets.  The  measured  EOD  coefficient  slightly 
increases  from  0.47  to  0.48  by  increasing  the  stack  current  from  1 
to  3  A,  which  is  in  accordance  with  the  data  of  Husar  et  al.  [  1 1  ],  but 
their  test  was  performed  at  higher  current  densities  and  using  a 
thicker  membrane. 

To  determine  a  specific  heat  capacity  for  the  stack,  a  constant 
current  is  drawn  for  a  short  period  of  time,  and  the  stack  temper¬ 
ature  evolution  is  measured.  The  test  resulted  in  a  specific  heat 
capacity  of  1260Jkg-1  K-1.  A  similar  value  is  also  used  in  the  heat 
transfer  model  of  He  et  al.  [18]. 

5.  Results  and  discussion 

5.1.  Steady-state  validation 

By  validating  the  steady-state  model,  the  general  model  settings 
such  as  the  boundary  conditions,  applied  physical  phenomena  or 
experimentally  determined  coefficients  can  be  checked.  Therefore, 
experimental  data  of  the  H-100  fuel  cell  stack,  directly  obtained 
from  laboratory  tests,  were  compared  to  the  simulation  results. 
The  initial  conditions  were  a  humidified  anode  and  cathode,  and 
a  constant  fan  flow  rate.  Even  though  the  regular  working  condi¬ 
tions  of  the  H-100  suggest  a  dry  anode,  these  tests  were  performed 
in  order  to  check  the  model  behavior  even  under  conditions  that  do 
not  appear  in  a  normal  operation.  The  stack  current  was  stepwise 
increased  from  1  to  3  A.  Fig.  13  shows  the  measured  and  modeled 
values  of  the  cathode  and  anode  outlet  water  concentrations,  as 
well  as  the  fuel  cell  stack  temperature  at  different  currents. 

5.2.  Dynamic  validation 

The  model  has  also  been  compared  to  dynamic  experiments  per¬ 
formed  with  the  studied  PEM  fuel  cell  stack.  These  tests  have  shown 
that  the  model  is  able  to  give  a  proper  dynamic  representation  of 
the  actual  stack  behavior.  Fig.  14  shows  the  comparison  of  the  sim¬ 
ulated  cathode  outlet  water  concentration  with  the  measurement 
data,  when  changing  the  cooling  fan  velocity  at  a  constant  working 
point  of  5  A. 


The  test  results  show  that  there  is  a  delay  in  the  cathode  outlet 
dew  point  measurement,  which  is  due  to  the  response  time  of  the 
dew  point  sensor  used  in  this  experiment.  However,  disregarding 
this  unavoidable  measurement  error,  the  simulated  cathode  outlet 
water  concentration  gives  a  good  representation  of  the  experimen¬ 
tal  results. 

The  second  validation  parameter  is  the  stack  temperature, 
which  was  also  measured  and  simulated  within  the  same  test.  The 
comparison  is  illustrated  in  Fig.  15. 

A  delay  in  the  measured  stack  temperature  can  be  noticed  here 
as  well.  This  is  not  only  due  to  the  response  time  of  the  sensor 
but  also  due  to  the  differences  between  the  programmed  velocity 
curve  and  the  measured  curve,  which  shows  a  smaller  slope  at  low 
velocities  and  therefore  the  temperature  changes  slower  than  in  the 
model.  Regarding  the  temperature  gradients,  the  model  matches 
the  experiment,  especially  during  the  cooling  phase.  The  modeled 
fuel  cell  stack  temperature  increases  faster  than  the  measured  tem¬ 
perature  when  the  fan  velocity  is  reduced.  The  slower  increase  in 
the  measured  temperature  could  be  caused  by  small  amounts  of 
convective  and  radiative  heat  removal  from  the  stack,  which  is  not 
included  in  the  model. 

The  conclusion  of  both  tests  is  that  the  model  demonstrates  an 
accurate  dynamic  representation  of  the  fan. 

5.3.  Sensitivity  analysis 

In  order  to  quantify  the  sensitivity  of  the  input  variables  of  the 
model,  a  sensitivity  analysis  was  performed  at  different  current 
set  points.  Table  5  shows  the  model  input  variable  settings  and 


Measured  inlet  air  velocity 

—  Programmed  inlet  air  velocity 

—  Measured  stack  temperature 

Simulated  stack  temperature 

Fig.  15.  Dynamic  validation— temperature  fan  test. 
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Table  5 

Input  and  output  variable  values  and  input  variable  variations. 


Tea  (°C) 

Tan  (°C) 

*H20,ca 

XH20,an 

vca  (ms-1) 

htack  (A) 

Vstack  (V) 

in 

30.00 

30.00 

0.01475 

0.23937 

0.4 

3.0 

13.5 

out 

35.84 

38.10 

0.01577 

0.50157 

- 

- 

- 

A$in 

1.0 

- 

0.00096 

0.01424 

0.1 

0.1 

0.1 

Table  6 

Sensitivity  matrix  at  3  A. 

tin 

Tca,out  (°C) 

S(Tca,out)  (%) 

Tan, out  (°C) 

s(Tan>0Ut)  (%) 

<^H20  ,ca,out 

s(Xh20 ,ca,out)  (Z>) 

Xh20, an,  out 

s(Xh20, an, out)  ('O 

Tca,in 

36.86 

2.85 

39.10 

2.62 

0.015788 

0.13 

0.502341 

0.15 

^H20  ,ca,in 

35.84 

0.00 

38.10 

0.00 

0.016744 

6.20 

0.509306 

1.54 

*H20  ,an,in 

35.84 

0.00 

38.10 

0.00 

0.015789 

0.14 

0.509386 

1.56 

Vca, in 

34.84 

-2.80 

37.30 

-2.10 

0.015570 

-1.25 

0.499355 

-0.44 

htack 

36.04 

0.56 

38.37 

0.71 

0.015802 

0.22 

0.510875 

1.85 

Vstack 

35.79 

-0.14 

38.03 

-0.18 

0.015765 

-0.02 

0.501492 

-0.02 

the  resulting  output  variable  values  at  a  stack  current  of  3  A.  It  also 
includes  the  variable  variations  for  the  sensitivity  analysis.  The  vari¬ 
ation  of  the  anode  and  cathode  inlet  water  mass  fraction  represents 
a  variation  in  the  measured  inlet  dew  point  temperature  of  1  °C, 
from  20  to  21  °C  at  the  cathode  and  from  25  to  26  °C  at  the  anode. 
The  big  difference  between  anode  and  cathode  water  mass  frac¬ 
tion  is  not  only  due  to  the  difference  in  dew  point  temperature  but 
rather  due  to  the  much  higher  gas  flow  rate  at  the  cathode.  Since 
the  anode  inlet  temperature  has  no  effect  it  is  not  included  in  this 
analysis. 

By  changing  the  input  variables  (§,•„),  according  to  the  variable 
variations  (zAf  defined  in  Table  5,  the  impact  of  input  signal  off¬ 
sets  or  variable  deviations  on  each  output  variable  can  be  observed. 
Table  6  shows  the  sensitivity  matrix,  depicting  the  resulting  output 
variable  values  for  each  input  variable  variation  with  the  respective 
relative  sensitivity  s. 

As  a  result  it  can  be  seen  that  increasing  the  cathode  inlet  dew 
point  temperature  by  1  °C  has  almost  the  same  effect  on  the  anode 
outlet  water  mass  fraction  as  increasing  the  anode  inlet  dew  point 
temperature,  however  not  vice  versa.  Also  the  stack  current  has  a 
similar  effect  on  the  anode  outlet  water  mass  fraction.  Moreover, 
an  increase  of  the  cathode  inlet  air  velocity  of  0.1  ms-1  has  the 
same  effect  as  a  1  °C  change  in  the  ambient  temperature  at  these 
conditions.  These  observations  indicate  that  the  cathode  mass  flux 
dominates  the  energy  and  mass  balance  in  this  stack.  All  variables 
show  very  low  sensitivity  to  a  variation  of  the  stack  voltage,  but 
are  sensitive  to  a  change  in  the  stack  current.  Repeating  the  exper¬ 
iment  at  higher  and  lower  stack  currents  leads  to  similar  variable 
sensitivities. 

6.  Conclusion 

A  two-dimensional,  non-isothermal,  dynamic  model  of  a  100  W 
open  cathode,  self-humidified  PEM  fuel  cell  system  has  been 
developed  with  respect  to  water  and  heat  transport  within  the 
cell.  The  crucial  coefficients  for  water  transport,  namely  the 
diffusion  and  the  EOD  coefficient,  have  been  determined  exper¬ 
imentally.  The  diffusivity  of  water  vapor  through  the  membrane 
at  30  °C  was  determined  to  be  3.3  x  10-8  m2  s-1  and  increases  by 
3  x  10-10  m2  s-1  °C_1  with  increasing  temperature  to  50  °C.  The 
EOD  coefficient  was  found  to  be  0.47-0.48  water  molecules  per  pro¬ 
ton  at  stack  currents  from  1  to  3  A.  Moreover,  the  bulk  specific  heat 
capacity  for  one  cell  unit,  consisting  of  a  MEA  plus  a  bipolar  plate, 
was  experimentally  determined  to  be  1260Jkg_1  K-1.  The  model 
has  been  validated  by  using  experimental  data  directly  obtained 
from  laboratory  tests  with  the  investigated  fuel  cell  stack,  which 


has  shown  that  the  model  predictions  match  the  experimental  data 
well.  The  model  is  kept  simple  and  is  capable  of  representing  system 
specific  control  mechanisms  for  water  and  heat  management,  as 
demonstrated  within  the  dynamic  validation.  As  it  combines  most 
of  the  physical  phenomena  that  occur  within  a  PEM  fuel  cell,  it 
permits  for  a  comprehensive  study  of  these  control  mechanisms. 
However,  the  model  can  still  be  improved  by  including  charge 
transfer,  two  phase  flow  characteristics  as  well  as  temperature 
driven  water  transport.  Moreover,  further  experiments  are  needed 
to  observe  the  dynamic  effect  of  water  storage  in  the  membrane  and 
the  anode  GDL,  and  finally  include  them  in  the  model.  The  devel¬ 
oped  model  is  intended  to  be  used  to  simulate  and  study  the  effects 
of  water  transport  and  its  influence  on  the  system  performance, 
and  to  develop  new  water  management  control  strategies,  that  are 
strongly  demanded,  as  recent  papers  have  shown.  The  model  is  easy 
to  handle  by  the  user-friendly  CFD  software  COMSOL  Multiphysics, 
and  can  be  easily  extended.  Furthermore,  it  is  applicable  to  other 
PEM  fuel  cell  systems,  following  the  developed  modeling  strategy 
and  performing  the  experiments  in  order  to  determine  the  specific 
coefficients. 
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